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Over five years we have compared the hyperfine frequencies of ^'^''Cs and *^Rb atoms in their 
electronic ground state using several laser cooled ^^'^Cs and ^'^Rb atomic fountains with an accuracy 
of ~ 10~^^. These measurements set a stringent upper bound to a possible fractional time variation 

d 1„ ( t^Rb 



of the ratio between the two frequencies : ^ In f ^ ) = (0.2 ± 7.0) x 10 



"i^yr ^ (la uncertainty). 

The same limit applies to a possible variation of the quantity (/iRb/MCs)^"'''^'', which involves the 
ratio of nuclear magnetic moments and the fine structure constant. 
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Since Dirac's 1937 formulation of his large number hy- 
pothesis aiming at tying together the fundamental con- 
stants of physics jl] , large amount of work has been de- 
voted to test if these constants were indeed constant over 
time |,|. 

In General Relativity and in all metric theories of grav- 
itation, variations with time and space of non gravita- 
tional fundamental constants such as the fine structure 
constant a = e"^ /Aireohc are forbidden. They would vi- 
olate Einstein's Equivalence Principle (EEP). EEP im- 
poses the Local Position Invariance stating that in a lo- 
cal freely falling reference frame, the result of any lo- 
cal non gravitational experiment is independent of where 
and when it is performed. On the other hand, almost all 
modern theories aiming at unifying gravitation with the 
three other fundamental interactions predict violation of 
EEP at levels which are within reach of near-future ex- 
periments [0, ||. As the internal energies of atoms or 
molecules depend on electromagnetic, as well as strong 
and weak interactions, comparing the frequency of elec- 
tronic transitions, fine structure transitions and hyperfine 
transitions as a function of time or gravitational potential 
provides an interesting test of the validity of EEP. 

To date, very stringent tests exist on geological and 
cosmological timescales. The analysis of the Oklo nu- 
clear reactor showed that, 2 x 10^ years ago, a did not 
differ from the present value by more than 10^^ of its 
value 1^. Light emitted by distant quasars has been 
used to perform absorption spectroscopy of interstellar 
clouds. For instance, measurements of the wavelengths of 
molecular hydrogen transitions test a possible variation 
of the electron to proton mass ratio me/rrip Q. Compar- 
isons between the gross structure and the fine structure 
of neutral atoms and ions would indicate that a for a 
redshift z 1.5 (~ 10 Gyr) differed from the present 
value: Aa/a = (-7.2 ±1.8) x 10"*^ @. Today this is the 
only claim that fundamental constants might change. 



On much shorter timescales, several tests using fre- 
quency standards have been performed [|| |l^, |ll| . These 
laboratory tests have a very high sensitivity to changes in 
fundamental constants. They are repeatable, systematic 
errors can be tracked as experimental conditions can be 
changed. 

In this letter we present results that place a new strin- 
gent limit to the time variation of fundamental constants. 
By comparing the hyperfine energies of ^'^^Cs and ®^Rb in 
their electronic ground state over a period of nearly five 
years, we place an upper limit to the rate of change of 
the ratio of the hyperfine frequencies v-Rh/vca- Our mea- 
surements take advantage of the high accuracy (~ 10~^^) 
of several laser cooled Cs and Rb atomic fountains. Ac- 
cording to recent atomic structure calculations |jll|, p^ , 
these measurements are sensitive to a possible variation 
of the quantity (^Rb/MCs)Q!~° '''*, where /i's are the nu- 
clear magnetic moments. We anticipate major advances 
in these tests using frequency standards, thanks to recent 
advances in optical frequency metrology using femtosec- 
ond lasers 

In our experiments, three atomic fountains are com- 
pared to each other, using a hydrogen maser (H-maser) 
as a flywheel oscillator (Figj^). Two fountains, a trans- 
portable fountain FOM, and FOl ||l^ are using cesium 
atoms. The third fountain is a dual fountain (DF) [ p^ , 
operating alternately with rubidium (DFRb) and cesium 
(DFcs)- These fountains have been continuously up- 
graded in order to improve their accuracy from 2 x 10~^^ 
in 1998 to 8 X 10"^^ for cesium and from 1.3 x lO^^'* |l|] 
to 6 X 10^^^ for rubidium. 

Fountain clocks operate as follows. First, atoms are 
collected and laser cooled in an optical molasses or in a 
magneto-optical trap in 0.3 to 0.6 s. Atoms are launched 
upwards, and selected in the clock level (rnp = 0) by a 
combination of microwave and laser pulses. Then, atoms 
interact twice with a microwave field tuned near the 
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hyperfine frequency, in a Ramsey interrogation scheme. 
The microwave field is synthesized from a low phase noise 
100 MHz signal from a quartz oscillator, which is phase 
locked to the reference signal of the H-maser (Fig.|l|). Af- 
ter the microwave interactions, the population of each 
hyperfine state is measured using light induced fluores- 
cence. This provides a measurement of the transition 
probability as a function of microwave detuning. Suc- 
cessive measurements are used to steer the average mi- 
crowave field to the frequency of the atomic resonance 
using a digital servo system. The output of the servo 
provides a direct measurement of the frequency differ- 
ence between the H-maser and the fountain clock. 
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FIG. 1: BNM-SYRTE clock ensemble. A single 100 MHz sig- 
nal from a H-maser is used for frequency comparisons and is 
distributed to each of the microwave synthesizers of the ^'^''Cs 
(FOl, FOM, DFcs) and *^Rb fountain clocks. In 2001, the 
Rb fountain has been upgraded and is now a dual fountain 
using alternately rubidium (DFab) or cesium atoms (DFcs). 
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FIG. 2: The 2002 frequency comparison data, a) H-maser 
fractional frequency offset versus FOM (■), and alternately 
versus DF^b (°) and DFcs (A between dotted lines), b) 
Fractional frequency differences. Between dotted lines, Cs- 
Cs comparisons, outside Rb-Cs comparisons. Error bars are 
purely statistical. They correspond to the Allan standard de- 
viation of the comparisons and do not include contributions 
from fluctuations of systematic shifts of Table ^ 



The three fountains have different geometries and op- 
erating conditions: the number of detected atoms ranges 
from 3 X 10^ to 2 x 10^ at a temperature of ^ 1 fiK, the 
fountain cycle duration from 1.1 to 1.6 s. The Ramsey 
resonance width is between 0.9 and 1.2 Hz. In measure- 
ments reported here the fractional frequency instability 
is (1 — 2) X 1G~^'^T~^/^, where r is the averaging time in 
seconds. Fountain comparisons have a typical resolution 
of ^ 10^^^ for a 12 hour integration, and each of the four 
data campaigns lasts from 1 to 2 months during which 
an accuracy evaluation of each fountain is performed. 

The 2002 measurements are presented in Fig.|[ which 
displays the maser fractional frequency offset, measured 
by the Cs fountains FOM and DFcs- Also shown is the 
H-maser frequency offset measured by the Rb fountain 
DFpib where the Rb hyperfine frequency is convention- 
ally chosen to be i/Rb(1999) = 6 834 682 610.904 333 Hz, 
our 1999 value. The data are corrected for the systematic 
frequency shifts listed in Table |[ The H-maser frequency 
exhibits fractional fluctuations on the order of 10"^'' over 
a few days, ten times larger than the typical statistical 
uncertainty resulting from the instability of the fountain 
clocks. In order to reject the H-maser frequency fluc- 
tuations, the fountain data are recorded simultaneously 
(within a few minutes). The fractional frequency dif- 
ferences plotted in Fig.|3b illustrate the efficiency of this 



rejection. DF is operated alternately with Rb and Cs, al- 
lowing both Rb-Cs comparisons and Cs-Cs comparisons 
(central part of Fig.0) to be performed. 



TABLE L Accuracy budget of the fountains involved in the 
2002 measurements (DF et FOM). 



Fountain 


DFcs 1 DFRb 1 FOM 


Effect 


Value & Uncertainty (10^^") 


2"" order Zeeman 


1773.0 ± 5.2 


3207.0 ± 4.7 


385.0 ± 2.9 


Blackbody Radiation 


-173.0 ±2.3 


-127.0 ± 2.1 


-186.0 ±2.5 


Cold collisions 
+ cavity pulling 


-95.0 ±4.6 


0.0 ± 1.0 


-24.0 ±4.8 


others 


0.0 ±3.0 


0.0 ±3.0 


0.0 ±3.7 


Total uncertainty 


8 


6 


8 



Systematic effects shifting the frequency of the foun- 
tain standards are listed in Table |. The quantization 
magnetic field in the interrogation region is determined 
with a 0.1 nT uncertainty by measuring the frequency of 
a linear field-dependent "Zeeman" transition. The tem- 
perature in the interrogation region is monitored with 
5 platinum resistors and the uncertainty on the black- 
body radiation frequency shift corresponds to tempera- 
ture fluctuations of about 1 K [p^. Clock frequencies 
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are corrected for the cold collision and cavity pulling fre- 
quency shifts using several methods All other 
effects do not contribute significantly and their uncer- 
tainties are added quadratically. We searched for the 
influence of synchronous perturbations by changing the 
timing sequence and the atom launch height. To search 
for possible microwave leakage, we changed the power 
(x9) in the interrogation microwave cavity. No shift was 
found at a resolution of 10~^^. The shift due to residual 
coherences and populations in neighboring Zeeman states 
is estimated to be less than 10~^^. As shown in the 
shift due to the microwave photon recoil is very similar 
for Cs and Rb and smaller than -1-1.4 x 10~^^. Relativis- 
tic corrections (gravitational redshift and second order 
Doppler effect) contribute to less than 10"^^ in the clock 
comparisons. 

For the Cs-Cs 2002 comparison, we find: 



(2002) 



,,FOM 



(2002) 



-12(6)(12) X 10" 



-16 



(1) 



where the first parenthesis reflects the la statistical un- 
certainty, and the second the systematic uncertainty, ob- 
tained by adding quadratically the inaccuracies of the 
two Cs clocks (see Table ^. The two Cs fountains are in 
good agreement despite their significantly different oper- 
ating conditions (see Table |), showing that systematic 
effects are well understood at the 10~^^ level. 

In 2002, the ^^Rb frequency measured with respect to 
the average ^'^^Cs frequency is found to be: 

i'Rb(2002) = 683468261G.904324(4)(7)Hz (2) 

where the error bars now include DFRb , DFcs and FOM 
uncertainties. This is the most accurate frequency mea- 
surement to date. 

In Fig.^ are plotted all our Rb-Cs frequency compar- 
isons. Except for the less precise 1998 data [|l^, two Cs 
fountains were used together to perform the Rb measure- 
ments. The uncertainties for the 1999 and 2000 measure- 
ments were 2.7 x 10^^*^, because of lower clock accuracy 
and lack of rigorous simultaneity in the earlier frequency 
comparisons |16 . A weighted linear fit to the data in 
Fig.^ determines how our measurements constrain a pos- 
sible time variation of t'Rb/'^Cs. We find: 



dt 



^ln( ^) =(0.2 ±7.0) X 10-16 yr"' 



(3) 



which represents a 5-fold improvement over our previous 
results and a 100-fold improvement over the Hg+-H 
hyperfine energy comparison pl . 

We now examine how this result constrains possible 
variations of fundamental constants. For an alkali with 
atom number Z , the hyperfine transition frequency can 
be approximated by: 



(4) 



MJD 

50500 51000 51500 52000 52500 53000 



(U 
(L> 




1997 1998 1999 2000 2001 2002 2003 2004 

Year 



FIG. 3: Measured *^Rb frequencies referenced to the ^'^'^Cs 
fountains over 57 months. The 1999 measurement value 
(i^Rb(1999) = 6 834 682 610.904 333 Hz) is conventionally used 
as reference. A weighted linear fit to the data gives 

at V "^Cs J 

spond to the la slope uncertainty. 



(0.2 ± 7.0) X 10"^''yr"\ Dotted lines corre- 



where i?oo is the Rydberg constant, c the speed of light, 
fi the magnetic moment of the nucleus, fijy the nu- 
clear magneton. Frei{Za) is a relativistic function which 
strongly increases with Z |l^, For ^'^'^Cs, this Casimir 
relativistic contribution amounts to 40 % of the hyper- 
fine splitting and a 



dln{F,.ai{Za)) _ 
da 



0.74. For 87Rb, this 



quantity is 0.30 |2^. Following |11|] and neglecting possi- 
ble changes of the strong and weak interactions affecting 
/iRb and ^cs, the sensitivity of the ratio J^Rb/j^Cs to a 
variation of a is simply given by: 



d 



a In, 



In 



(0.30 - 0.74) 



-0.44. 



(5) 



Using equations ^ and |^, we thus set the new limit: 
a/a = (-0.4 ± 16) x lO-^^yr^i. 

In contrast with [|ll] , Ref . |2|] argues that a time varia- 
tion of the nuclear magnetic moments must also be con- 
sidered in a comparison between hyperfine frequencies. 
The magnetic moments /i can be calculated using the 
Schmidt model. For atoms with odd A and Z such as 
*^Rb and ^^'^Cs, the Schmidt magnetic moment /i'*' is 
found to depend only on gp, the proton gyromagnetic 
ratio. With this simple model, Ref. p2[ finds: 




(6) 



Attributing any variation of v-Rh/i'Cs to a variation of 
gp, equations || and ^ lead to: gp/gp = (0.1 ± 3.5) x 



10"^'' yr"^. However, it must be noted that the Schimdt 
model is over simplified and does not agree very accu- 
rately with the actual magnetic moment. 
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Moreover, attributing all the time variation of t'Rb/t'Cs 
to either gp or a independently is somewhat artificial. 
Theoretical models allowing for a variation of a also al- 
low for variations in the strength of the strong and elec- 
troweak interactions. For instance, Ref. j|] argues that 
Grand unification of the three interactions implies that 
a time variation of a necessarily comes with a time vari- 
ation of the coupling constants of the other interactions. 
Ref. j|] predicts that a fractional variation of a is ac- 
companied with a ~ 40 times larger fractional change 
of irie/mp. In order to independently test the stability of 
the three fundamental interactions, several comparisons 
between different atomic species and/or transitions are 
required. For instance and as illustrated in , absolute 
frequency measurements of an optical transition is sensi- 
to a different combination of fundamental constants : 
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(A*Cs /A<jv)('iT-e/™p)'^'^' where x depends on the particular 
atom and/or transition. 

A more complete theoretical analysis going beyond the 
Schmidt model would clearly be very useful to interpret 
frequency comparisons involving hyperfine transitions. 
This is especially important as most precise frequency 
measurements, both in the microwave and the optical 
domain ^ |2^, are currently referenced to the ^^'^Cs 
hyperfine splitting, the basis of the SI definition of the 
second. The H hyperfine splitting, which is calculable to 
a high accuracy, has already been considered as a possi- 
ble reference several decades ago. Unfortunately, despite 
numerous efforts, the H hyperfine splitting is currently 
measured to only 7 parts in lO^'^ (using H-masers), al- 
most three orders of magnitude worse than the results 
presented in this letter. 

In summary, by comparing ^'^^Cs and ^^Rb hyperfine 
energies, we have set a stringent upper limit to a possi- 
ble fractional variation of the quantity (/iRb//^Cs)Q:^*''^^ 
at (—0.2 ± 7.0) X 10~^^yr~^. In the near future, accura- 
cies near 1 part in 10^^ should be achievable in microwave 
atomic fountains, improving our present Rb-Cs compar- 
ison by one order of magnitude. 

We anticipate that comparisons between rapidly pro- 
gressing optical and microwave laser-cooled frequency 
standards currently developed in several laboratories will 
bring orders of magnitude gain in sensitivity. In order to 
have the full benefit of these advances, frequency com- 
parisons with improved accuracy between these distant 
clocks will be necessary. Serving this purpose, a new 
generation of time/frequency transfer at the 10~^^ level 
is currently under development for the ESA space mission 
ACES which will fly ultra-stable clocks on board the in- 
ternational space station in 2006 ||2^ . These comparisons 
will also allow for a search of a possible change of fun- 
damental constants induced by the annual modulation of 
the Sun gravitational potential due to the elliptical orbit 
of the Earth 
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